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We study gravitational collapse of a spherically symmetric inhomogeneous dust cloud in the Love- 
lock theory without cosmological constant. We show that the final fate of gravitational collapse in 
this theory depends on the spacetime dimensions. In odd dimensions the naked singularities formed 
are found to be massive. In the even dimensions, on the other hand, the naked singularities are 
found to be massless. We also show that the curvature strength of naked singularity is independent 
of the spacetime dimensions in odd dimensions. However, it depends on the spacetime dimensions 
in even dimension. 
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I. INTRODUCTION 



Whether gravitational collapse results in a black hole 
or a naked singularity is one of the crucial issues in clas- 
sical general relativity. A singularity is not a "point" 
in the spacetime and information from such a place is 
in fact an end of causality in the theory. In order to 
avoid such a breakdown of physics, it is often believed 
that naked singularities cannot be formed under physi- 
cally reasonable conditions in classical general relativity 
(GR). This is known as the cosmic censorship conjecture 
(CCC) Q. Broadly there are two versions of the CCC, 
i.e., a strong cosmic censorship conjecture (SCCC) and a 
weak cosmic censorship conjecture (WCCC). The SCCC 
prohibits the formation of locally naked singularity, i.e., 
singularities are not visible to even nearby observers. On 
the other hand, the WCCC prohibits the formation of 
globally naked singularity, allowing a local observer mov- 
ing around singularity to see it. Though there have been 
many efforts to prove the CCC, it remains as one of the 
important unsolved problem in GR. Meanwhile, a con- 
siderable number of counterexamples have been found 
violating both the WCCC and SCCC. 

In 1939 Oppenheimer and Snyder studied the gravita- 
tional collapse of a spherically symmetric homogeneous 
dust cloud and found that the collapse ends in a black 
hole 1^. However, in absence of a general proof of the 
CCC there have been studies on collapse using various 
matter models like dust, radiation, scalar fields, and so 
on. In these exact solutions of Einstein equations several 
counterexamples to the CCC have been found. For exam- 
ple, naked singularities form in the collapse of spherically 
symmetric inhomogeneous dust models [l-il (see Ref. 
p^ft] and the references therein) . In general there are two 
types of singularities in spherical collapse, namely, shell- 
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crossing and shell-focusing singularities. But, in what 
follows, we will focus only on the central shell-focusing 
singularities. 

Recently, higher-dimensional spacetimes have gener- 
ated considerable interest in theoretical physics. These 
higher-dimensional models are motivated by string/M 
theory. The case of collapse of a inhomogeneous dust 
cloud in the simple extension of GR to higher dimensions 
was well studied. It was shown there that the singularity 
cannot be naked when the spacetime dimension is more 
than six for smooth initial data [IMl (see also Ref. [l^ 
and the references therein). More precisely, one could 
confirm that the SCCC holds in dimensions higher than 
six in such a naive extension of GR to higher dimensions. 
Clearly, the spacetime dimensions play an important role 
in the final fate of the gravitational collapse. However, 
it is expected that spacetime should be extremely curved 
near singularities putting in doubt the validity of classical 
GR in these regions. Indeed the string/M theory predicts 
the higher curvature corrections to the Einstein equation 
in its low-energy limit. In view of the above facts it is im- 
portant to study gravitational collapse in such a theory. 
Though the CCC was originally formulated for general 
relativity, it is worthwhile to investigate the implications 
of stringy effects on it. In this paper, we shall address 
this issue while focusing on the gravitational collapse of 
spherical symmetric dust clouds. But there is no canon- 
ical recipe to include higher-order curvature corrections 
to GR. Although string theory predicts such corrections, 
we do not have a clear picture yet. In this paper we shall 
adopt a special combination of higher curvature correc- 
tions so that the field equations do not include the higher 
derivative terms of the metric than the third order, that 
is, the Lovelock theory [14:]. There is no strong reason 
why the CCC should hold in the Lovelock theory, since 
energy conditions are generally violated in this theory. 
Nevertheless, if singularities do form in this theory, it is 
important to know their features. 

In the Lovelock theory the static black hole solutions of 
a spherically symmetric vacuum were found in Refs. [H- 
[iTj . In Ref. [H, [l^ the stability of a static vacuum black 
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hole solution was discussed. The gravitational collapse of 
a spherical homogeneous dust cloud was analyzed in the 
dimensionally continued gravity [lO, [HI (the restricted 
class of the Lovelock theory [1^). It was shown there that 
the singularity will be covered by event horizon as in the 
Oppenheimer-Snyder's case (H. Subsequently, the grav- 
itational spherical collapse of inhomogeneous dust cloud 
in the Gauss-Bonnet theory of gravity was addressed in 
Ref. [l^l and it was shown that the final fate of collapse 
depends on the dimensions of spacetime. An exact so- 
lution and complete analjysis in the 5-D Gauss-Bonnet 
theory was given in Rcf. [24]. In five dimensions, singu- 
larities were found to be naked and massive. From six to 
nine dimensions, singularities are naked and massless. In 
dimensions higher than ten, singularities are always cov- 
ered. Therefore, including the Gauss-Bonnet corrections 
to GR makes the final state of collapse rather nontrivial. 
The gravitational collapse of null dust fiuid was also stud- 
ied [25-27j. These works also showed that the Lovelock 
corrections modify the collapse scenario in a nontrivial 
way. 

In this paper we investigate the gravitational collapse 
of a spherical inhomogeneous dust cloud in the Lovelock 
gravity, without cosmological constant and in any dimen- 
sions. We show that singularities do form in the final 
stage of gravitational collapse in general and that they 
could be naked for some initial dates. Comparing with 
the Gauss-Bonnet gravity, we found that the Lovelock 
term changes the nature of singularities. In the Gauss- 
Bonnet case, singularities which appear in more than six 
dimensions are massless. In the Lovelock theory with odd 
dimensions, we see that singularities are massive. On the 
other hand, singularities are massless in the even dimen- 
sions. This is an unique feature of the Lovelock gravity. 

The rest of this paper is organized as follows. In Sec. 
II, we briefly review the Lovelock theory. In Sec. Ill, 
we derive the basic equations for gravitational spherical 
dust collapse in the Lovelock theory. In Sec. IV, we ana- 
lyze the nature of the singularity and apparent horizon. 
In Sec. V, we show the existence of null geodesies from 
the singularities. In Sec. VI, we analyze the strength of 
the naked singularities. Finally we will give our conclu- 
sion and discussion in Sec. VII. In the Appendix A, we 
discussed the junction between inner and outer solutions. 
In the Appendix B, the analysis of the homogeneous case 
will be briefly given. In this paper, we employ the units 
of c = 1 and SttG = 1, where c is the speed of light and 
G is the gravitational constant in higher dimensions. 



II. LOVELOCK THEORY 

In this section we give an overview of the Lovelock 
theory m D = n-\- 2-dimensional spacetimes. This is the 
most general theory of gravity satisfying following three 
conditions: 

1. The field equations are written in terms of a sym- 
metric rank-2 tensor. 



2. The theory is consistent with the conservation law 
of the energy-momentum tensor. 

3. The theory does not include higher than third order 
derivatives. 

The Lagrangian of the theory is given by 



£ - V —C 
^-^ m 



m=l 



(1) 



where 



■ ^^'^2m-l'^2m 



(2) 

In the above, R^u'^'^ is the Riemann curvature tensor, 
and '5j^i^2'.'.'.'^^2'^"l\'/i^l is the generalized and totally anti- 
symmetric Kronecker delta. {a„i} are arbitrary constants 
which cannot be determined by the theory itself. The 
suffixes /ii . . . /i2m and vi . . . V2m run from 1 to Z), and 
/c is a constant depending on the spacetime dimensions, 
defined hy k — [{D — l)/2] {[x] is the integer part of x). 
Throughout this paper we suppose that a™ are positive. 
The field equations for this Lagrangian are 
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MMlA'2---A'2m-lA'2n 



Y R ^1/^2 D ^2711-1/^217 
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(3) 



where is energy-momentum tensor of matter. It 
should be remembered that V^tJ^ holds, where 
is the covariant derivative with respect to the metric g^^. 

There is an exact solution of the vacuum, static and 
spherical symmetric black hole [T5l - [T7| . The metric is 
given by 



1 



1^2 



dx^dx-' , 



(4) 



where 7^ is the line element of n-dimensional constant 
curvature hypersurface and /(r) is 



(5) 



In the above, k is the constant curvature of the n- 
dimensional hypersurfaces which take values —1,0 and 1. 
In this paper we are interested only in asymptotically flat, 
spherically symmetric spacetime and then set k = I. The 
function tp{r) is a solution to the algebraic equation 



m 



2m-2 

n (.^-p)\r 

p=l 



^n+1 ' 



(6) 



Here /i is a constant proportional to the ADM mass. 

Under the assumption that all the coefficient 
positive, W^['0] is a monotonically increasing function of V' 
when > 0. Therefore, positive V' has a unique solution 
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to Eq. ([6]). When we take r 
goes to zero as 



oo, the positive solution 



+ (7) 

From this fact and the metric form of Eq. ([5|), it is easy 
to see that the solution corresponds to an asymptotically 
flat spacetime. A negative V' could also be the solution. 
But in this case, the solution corresponds to an asymp- 
totically anti-deSitter (AdS) spacetime. In this paper we 
shall consider on the positive ip solution because we are 
interested in asymptotically flat spacetimes. 

The horizon is located at / = 0. Then the equation 
for the horizon is 



1 



and Eq. 1^ implies 



(8) 



(9) 



at the horizon. The solution to Eq. (jH) depends on space- 
time dimension, n. In even dimensions, Eq. Q becomes 



n/2 |- , 2m-2 

m=2 L ^ p=l 



(n-p) }tp' 



1 = 0. 



(10) 



It is easy to see that there are always solutions to Eq. 
(fTOl) . In odd dimensions, on the other hand, Eq. (|9]) 
becomes 



.2a, 



(«-l)/2 



p=l 

(n-l)/2 p , 2m-2 



m — 1 



-1-1 =0.(11) 



There are the solutions to the equation above if 
2a ,.4- '"-1 



(12) 



is satisfied. In odd dimensions, there is a tendency that 
spacetimes with small mass have no horizons. This is 
the characteristic feature which is different from even- 
dimensional cases. For example, in the five-dimensional 
asymptotically flat solution, f{r) becomes 



The the event horizon is given by 



r = - 02. 



(13) 



(14) 



Since 02 > 0, the size of the Gauss-Bonnet black hole is 
smaller than that of the corresponding five-dimensional 
Schwarzschild black hole. Therefore, the Gauss-Bonnet 
correction weaken the gravity. Note that this solution 
has the naked singularity when < 02. 

For our current study, the vacuum exact solution men- 
tioned above stands for the outer region of the collapsing 
dust cloud. Actually, we can join this outer solution with 
the inner cloud solution (see the appendix A for the de- 
tails). 



III. SPHERICAL COLLAPSE OF DUST 
CLOUDS 

In this section, we will consider the spherical collapse 
of dust clouds in the Lovelock theory. We first write down 
the basic equations. 

The energy-momentum tensor of dust cloud is 



(15) 



where s{t,r) is the energy density, which is supposed to 
be positive, and is the velocity. The metric is 

ds^ = -B{t, rfdf + A{t, rfdr^ + R{t, rfdill. (16) 

Here A{t, r), B{t, r) and i?(i, r) are arbitrary functions of 
t and r, and dfi^ is the line element of the n-dimensional 
unit sphere. Using the gauge freedom of coordinate r, 
we can always set i?(0,r) = r. In the current paper, we 
assume dtR < at initial surface. This means that the 
dust cloud is initially collapsing. 

From the conservation equation of the energy- 
momentum tensor 



v.t: = 0, 



we have two equations 



(17) 

(18) 
(19) 



In the above the dot and prime denote the partial deriva- 
tive of t and r, respectively. 

It is easy to see that Eq. implies B = B{t). Then, 
without a loss of generality, we can set B ^ \. Now the 
metric becomes 




ds^ 



-df + A{t, rfdr^ + R{t, rfdVL. 



2 jo2 

71 1 



(20) 



and the gravitational field equations for this simplified 
metric are 
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k 2m-2 

^ a™ n ~ ^) 



R 
R 



RA 



1 \ I R (n- (2m- 1)) / / R 
^ \r^ 



RAJ ^ i?2 



0,(21) 



k 2m-2 

n ('^-^') 1(9)- 

m— 1 p—0 



i?" R'A' 



m — 1 



RA 



1 

i?2 



i?A (n - (2m - 1)) / / i? 
RA ^ 2^^^^ vl ^ 



RA 



^1 -(^) +]^) 1=^(^.0, 



(22) 



and 



k 2m -2 



E 

m— 1 p—0 



n 



RA ^ i?2 



i?' R'A 



RA^ RA^ 



= 0. 



(23) 



From Eq. ([23]), we see that 

R' R'A 



= 



holds. This equation can be easily integrated as 
A ^' 



(24) 



where W{r) is the arbitrary function of r. Introducing L 
by 



2m -2 



2m 

m—1 p—0 



171 ^ 



i?2^ 



(26) 



and using Eq. (|25|) . we can rewrite Eqs. pip and ((22|) in 
the following simple forms 



J_d_ 

RR^ dt 
1 d 



[LR"+^) = 0, 
:ii?"+i) =e(i,r) 



Then the integration of the above equations implies 



_ n F{r) 

Lj - — 



2 ' 



(27) 
(28) 

(29) 



where F[r) is a function of r, which can be interpreted 
as a sort of mass function defined as 



F{r) ^^-f 
n Jo 



eR^'R'dr. 



(30) 



Indeed, F(r) is proportional to the generalized Misner- 
Sharp mass [2^1 ■ Using F{r), Eq. ()28p can be reexpressed 
as 



n F(r)' 
2R^R'' 



(31) 



For later convenience, we shall introduce the coeffi- 
cients ci by 



1 



if ; = 1 



ltnpLt("-p) if 2<;<fc. 



(32) 



(25) Using this, Eq. (f29|) becomes 



^ V / i?2'n Rn+1 

m—1 



(33) 



This is the key equation for the gravitational collapse. 

From Eq. (ISTI) . we can see two types of singularities 
which correspond to blowing up of energy density, that is, 
the shell- focusing singularities which form at at i? = 0, 
and the shell-crossing singularities at i?' = 0. In this 
paper we assume i?' > to exclude the shell-crossing 
singularities from our consideration. 

To make the initial surface regular at the center (r = 
0), the function F{r) should behave as 



F{t) 



M{r) 



(34) 



where M(r) is some regular function on initial surface. 

For later convenience, we introduce the function v(t, r) 
through 



i?(t, r) = rv{t, r). 



(35) 



Then our current initial conditions of R{t, r) are ex- 
pressed as 



«(0,r) = 1, 
w(0,r) < 0. 



(36) 
(37) 



The singularities at i? = are located at the orbit satis- 
fying v(t, r) = 0. 
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Before closing this section, we have one remark. At 
first glance, 



m=l p—0 



(38) 

is also the solution to Eq. (P^ . Combining above with 
Eq. (HI]), we have 



k 2m-2 



a,n n i^-P) 

m—l p—0 
■ \ 2 



{n- {2m- 1)) 
2to 



Together with Eq. (|22|) . we can show that this implies 
e{t,r) = 0. This means that the spacetimes must be 
vacuum, which is a solution but not our interest here. 



spectively. For the current metric form, dR/dt\± becomes 

• dr 

= i? + 

+ 

= R + 



dR 

m 

dR 



dt 
dr 
lit 



R' = R+1, 
R' = R- 1, 



(42) 
(43) 



where we used the fact of dr/dt\± = ±R'~ . The ap- 
parent horizon is defined to be the boundary of trapped 
regions. In our case it is easy to see that the location of 
the apparent horizons can be found through the solving 
of 

g'^'Rf.R. = 0. (44) 
For the metric of Eq. ([20]) . this condition becomes 



R-" = 1. 



A. Even dimensions 



(45) 



IV. APPARENT HORIZONS AND 
SINGULARITIES 

In this section we analyze the structure of spacetimes 
obtained in the previous section. We shall focus on the 
behavior of trapped surfaces and apparent horizons sur- 
rounding the singularities. 

We shall investigate the even and odd-dimensional 
cases separately since it is clear from Eqs. ^TU\\ and PT|) 
that they should have quite different features. As stated 
in Sec. II, we assume that all coefficients c; are positive. 
Here we also take W{r) = 1, which corresponds to the 
"marginally bound case" . 

First of all, it should be noted that bounce {v — 
and V ^ 0) cannot occur in this case for initially collaps- 
ing configurations. Then we will be able to expect the 
appearance of the shell-focusing singularity in general. 
Using Eqs. §^ and Eq. (gS]) becomes 

^ /^V" Mir) 

E^"(-j (40) 

m— 1 ^ ^ 

Substituting v = into Eq. (|40l) . then we have 

Af (r) = (41) 

for the bounce. This contradicts our assumption that e 
is positive. From this fact, bounce cannot occur and the 
cloud continues to collapse. We can also see from this 
that the range of v is from to 1. 

A trapped surface is described by the surfaces which 
satisfy dR/dt\± < 0, where d/dt\^ and d/dt\- are deriva- 
tives along the outgoing and ingoing null geodesies re- 



We first consider the apparent horizon in even dimen- 
sions. Using Eq. dJS]), Eq. ([55]) implies 

k=n/2 

ffe.en(i?)= ^ C„i?"-2™+l = ^^(^) ^ (45) 
ni—1 

which determines the location of the apparent horizons. 
Note that the left-hand side of the equation above is an 
odd function of R, and F(r) is a monotonically increas- 
ing function of r, with limr_>.o F{r) = 0. We depict the 
behavior of the left-hand side of Eq. in Fig. [1] 

As shown in the figure, for any given finite r, we always 





/ F(r) 








R 



FIG. 1: The profile of function Heven{R)- The thick line 
is Heven{R) and the dotted line is F{r). Heven{R) is the 
odd function and we assume that all coefficients are positive. 
Therefore, Heven is monotonically increasing function for R > 
and Jye„en(0) = 0. 
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have a solution to Eq. p6|) . This imphes that the ap- 
parent horizons always exist in these spacetimes and all 
singularities which occur at finite r > (but, R ~ 0) 
will be the inside of the apparent horizon. Therefore, 
only the central singularity at r = could be naked. In 
addition, they are massless because F{0) — 0, in even di- 
mensions. (Note that F{r) is proportional to a quasilocal 
mass.) This is the same property as in the Gauss-Bonnet 
case. In the Gauss-Bonnet case, all singularities except 
for r = were found to be covered by apparent horizons 
in any dimensions higher than six [23l |. 

Next, let us confirm that the region inside the ap- 
parent horizon is trapped. From Eq. p3|) . contin- 
ues to increase during the collapse, which means R de- 
creases. This implies R is monotonically decreasing be- 
cause we assume that the dust cloud is initially collaps- 
ing. Then, we can see that the region inside of the appar- 
ent horizon(_R = —1) is actually trapped, dR/dt\± < 0, 
from Eqs. (02]) and (|43| . 



Odd dimensions 



As in even dimensional cases, Eqs. (|45p and ([33l) give 



fc=(n+l)/2 
Hodd{R) = ^ 



= F{r) , 



(47) 



m— 1 



which determines the location of the apparent horizons. 
Contrary to what we had in even dimensions, the left- 
hand side of Eq. (|47| is the even function of R (See Fig. 
[5]). In the range of < r < rah, we see that there are no 





/ F(r) 








R 



FIG. 2: The profile of the function Hadd(R)- The thick line is 
Hodd{R) and the dotted line is F{r). Hodd{R) is the even func- 
tion and we assume that all coefficients are positive. There- 
fore Hodd is monotonically increasing function of i? > and 
Hodd{0) = C(n+i)/2 > 0. The figure shows that r satisfying 
F{r) < C(„^i)/2 have no solution to Eq. (|47p . 



solutions, where rah is determined through 

F{rah) = C(„+i)/2. 



(48) 



Hence the singularities which occurs in the range < 
f < fah are not wrapped by the apparent horizons. They 
could be naked and massive ( F{r*) > where < r* < 
rah) in odd dimensions. This result is in agreement with 
what is observed in the five-dimensional Gauss-Bonnet 
case ^23,]. For the Gauss-Bonnet case in any dimensions 
higher than six, however, the singularities at finite r shall 
always be safely wrapped. We also see that the region 
inside of the apparent horizon is trapped region from the 
same argument in even-dimensional cases. 



V. NULL GEODESICS FROM SINGULARITIES 

In this section, we shall analyze whether singularities, 
which are not wrapped by apparent horizons, can be 
naked or not. To this end we examine the existence of 
future-directed outgoing radial null geodesies from these 
singularities. Again, we need to analyze odd- and even- 
dimensional cases separately. 

From the metric of Eq. (PU)) . the outgoing radial null 
geodesic equation becomes 



^ = R'. 

dr 



(49) 



To show the existence of solution for this differential 
equation from the singularity, we shall employ the fixed- 
point method [1, 0, [2J| . We define a new function 



R 



(r-r,)9' 



(50) 



where g is a positive constant and is the location of 
singularities. Using x, the null geodesic equation can be 
rewritten as 



dx qx 
dr r ^ rs 

where 5(a;, r) is given by 



S(a;,r) 



R' 



(r 



^{x,r), 



1 + R 



(51) 



(52) 



One may consider the cases in which the function S(a;, r) 
can be expanded around r = r^ 



go(g)a^ + E»^ig'(g)'^° 
r - rs 
+ oiir~r,r), 



E{x,r) 



hiq,x){r-r,r 
(53) 



where a (> —1) and ^ 1 are constants. go,gi and h 
are some functions. > 1 is an integer. We will be able 
to use the following useful lemma to show the existence 
of radial null geodesic from singularities 
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Lemma [Lemma 10 in Ref. f^j //- // S(a;, r) is ex- 
panded as Eq. i5S\) with N — 1 and goi<]) ^ q, there 
exists an asymptotic solution satisfying x{rs) = [51/(9 ~ 
go)]^^^^~"'''^ near r — rs, and moreover it is the unique 
solution to Eq. i51\) which is continuous atr — rs- 



By virtue of this lemma, it is enough to investigate the 
features of the function 2(a;,r) for our current purpose. 
Hereafter we are interested in the future-directed outgo- 
ing null geodesies, which correspond to positive ^(r^), as 
can be seen from Eq. (H^ . 

From now on, we will address the existence of the null 
geodesies in odd- and even-dimensional cases separately. 



A. Odd dimensions 

As seen in the previous section, odd-dimensional space- 
times may have complicated structures. 



1. Singularities at rs = Q 

We first analyze the geodesies from singularity at r = 
in odd dimensions. In odd dimensions, the gravitational 
Eq. (Uni) becomes 



(n+l)/2 ... 2m , w . 

El v\ M ir) 

m— 1 



V 



n+l 



(54) 



Solving this with respect to i;, we have the following for- 
mal solution 

v^Vo{v,r). (55) 
The formal integration of this with respect to v is 

ds 



t{v,r) 



1 Vois,r)- 



(56) 



Let us assume the smooth distribution of matter near 
the center(r = 0) as 



M{r) = Mo + A'hr^ + AUr'^ + . . 
Then t{v, r) behaves near r = as 



t{v,r)=t{v,0) + —X2{v) + 



(57) 



(58) 



The orbit of the singularities is given by ts{r) = t{0,r). 
Since the definition of X2{v) is complicated, we will not 
write down its explicit form. Regardless of the theory, it 
is shown that X2(0) is non- negative under the condition 
of i?' > and i? < in Ref. [2^ . Here we further assume 
that X2(0) is positive, that is, X2(0) cannot achieve zero 



(this corresponds to the assumption for the initial data) 
Using the relation, which always holds. 



we have 



v' = {-rx2{v) + . . .)v. 



(59) 



(60) 



Now we set g = 3 in Eq. (I52|) and then we can expand 
E{x, r), with the help of Eqs. ([Ml), ^ and as 



^(x,r) 



V + rv' 



x + X2iv) 



(1 + rv) 



Mo 



l/(n+l) 



+ 0{r°). (61) 

From this expansion, the existence of solution to the null 
geodesic Eq. (|5T]) with 



x(0) 



X2(0) / Mo 



l/(n+l) 



2 \ c. 



> (62) 



is guaranteed according to the lemma. Since the posi- 
tivity of x{0) means the null geodesic is future-directed 
and outgoing null, this shows that the singularities are 
at least locally naked. 



2. Singularities at < rs < rah 

Next let us consider the geodesies from singularities at 
< Ts < rah- We consider the initial condition of 



M(r) = Mo + Mi(r - r,) + M2{r - r,f + 



(63) 



In the previous case of = 0, we assumed smoothness 
at the center which implies no kinklike distribution of 
matter. Therefore the next-to-leading term of expansion 
is the order of 0{r^). However, r — rg is not the 
symmetric center and then the next-to-leading term of 
expansion will be 0(r — r^) in general. 

Near the singularities, then, t{v,r) behaves as 



t{v, r) = t{v, 0) + {r - rs)xi {v) + 



(64) 



Following Ref. [2J], it is shown that xi(0) is non- 
negative. Here we simply assume that xi(0) cannot 
achieve zero. Using Eqs. ((5^)) and ([M)) . we see 



v' = (-Xi(") + ■••)"■ 



(65) 



Setting q = 1 in Eq. (15^ . we can expand 5(i,r), with 
the help of Eqs. ([5^ 1) . (P]) and as 
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5(i,r) = 



-Xi(w) 



Mn 



l/(n+l) 



l-r. 



Mn 



l/(n+l)N 



+ O((r-r,)0). 



(66) 



According to Eqs. p8|) and (|63| . rah is 

-c„+i \ i/("+i) 



Mo 



(67) 



From this expansion, the existence of a solution to the 
null geodesic Eq. ([51]) from singularities at r = (0 < 
rs < rah) with 



Mo 



l/(n+l) 



l-r. 



Mo 



Cn + l 



l/(n+l)^ 



> 



(68) 



r 



is proved according to the lemma. This means that the 
singularities are at least locally naked. 



B. Even dimensions 

Finally we shall consider the even-dimensional space- 
times. In even dimensions the gravitation field equation 
of (|40|) becomes 



^ fvY"" _ M{r) 

m— 1 ^ ^ 

Then we obtain formally 

V = V;(w,r). 
The integration of this with respect to v is 

r ds 
/i K(s,r)- 



(69) 

(70) 
(71) 



We also consider the smooth distribution of matter as 
M(r) = Mo + Mir^ + M^r^ + . . . . (72) 
Then t(v^ r) behaves like 



i(w,r) =t(v,0)-Fx2(«)y + . 



(73) 



As in the case of odd dimensions, we assume that X2(0) 
is positive. According to Eq. ([59]) . we have 

v' ^{~rx2{v) + ■■■)v. (74) 

Setting q — ^^qrp in Eq. (15^ . we can expand S(a;,r), 
with the help of Eqs. ^ and dH]), as 

l/n 



5-f X2(t')S"'/" 
S(x,r) = + 0(r-2/("+i)) 



(75) 



From this expansion, the existence of solution to the null 
geodesic Eq. (ISTI) with 



„ , 1 \ "/(n + l) / \ 



2n 



X2(0)"/^"+^^ (^-^j > 

(76) 

is guaranteed according to the lemma. Then the singu- 
larity is at least locally naked in this case too. 



VI. STRENGTH OF SINGULARITY 

We finally compute the strength of naked singularity, 
which is the most important feature characterizing sin- 
gularities. We define 



7^ — Rui^k^ k , 



(77) 



where /c^ — is the tangential vector of future out- 
going null geodesies from singularity and A is the affine 
parameter of the geodesies. We will check the behav- 
ior of TZ near the singularity (See Ref. [l^ for a more 
systematic definition of the strength of singularity). 

Consider the future outgoing radial null geodesies from 
singularities. We have a relation from the metric of Eq. 
(pg| as follows 

fc* = (78) 
By a straightforward calculation using Eq. (1751) , we have 
nF{ry 



The radial null geodesic equation is 
dk"- fR" 



R'ik'-f 



dX 



R' 



2R' (k 



u'■^2 



0. 



(79) 



(80) 



We need to consider odd- and even-dimensional cases 
separately. 
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A. Odd dimensions 



1. Singularity at r — 



Near the singularity at r = 0, we can compute TZ as 

/ \2/(n+l) 

TZ ^3nif^\ {¥f. (81) 

The geodesic equation near the singularity becomes 



d\ r 

The integration of Eq. gives us 



(82) 



(83) 



We now evaluate the behavior of TZ near the singularities 
as 



7^- A~3. 



(84) 



In odd-dimensional cases, we confirmed that TZ diverges 
at r = and the behavior of divergence does not depend 
on the spacetime dimensions. 



2. Singularity at < rs < rah 
Near the singularities at r = rs, TZ behaves as 

2/(n+l) / 



„ nr, / Mn , 
TZ — — - 1 



n + 1 Mo 



The geodesic equation becomes 



dX 



+ C(F)"~0, 



(85) 



(86) 



where C is a constant depending on the initial data. 
From Eq. ((55)) we see that 

r~rs^X, (87) 

near the singularity (A ~ 0). Then we can see 

7^-A-^ (88) 

As a result, TZ diverges at r = Ts (0 < rg < rah) and the 
behavior of divergence does not depend on the spacetime 
dimensions. 

B. Even dimensions 

1. Singularity at r — 
We can evaluate TZ near the central singularity as 



7^- (3n+ 1) 



Mo 
x(0)c^ 



2/r: 



r^2 



(89) 



The geodesic equation near the singularity is 
dk'' 2n 



+ 



dX {n + l)r 
The integration of Eq. (IMl) implies 



{ky~o. 



r ^ A3-+1 . 



Thus, near singularity, we have 

TZ^X-<^^) 



(90) 



(91) 



(92) 



In the even-dimensional cases, TZ diverges at r = and 
its behavior depends on the dimensions. This feature is 
different from odd-dimensional cases. 



VII. SUMMARY AND DISCUSSION 

In this paper we considered gravitational collapse of a 
spherical inhomogeneous dust cloud in the Lovelock grav- 
ity with any dimension. We found that the formation of 
apparent horizon depends on the dimensions, that is, if it 
is odd or even. In even dimensions, noncentral singulari- 
ties (r 7^ 0) will be safely wrapped by apparent horizons 
and only "central" singularities at r = could be naked. 
In odd dimensions, on the other hand, even those singu- 
larities which form at < r < rah may not be wrapped 
by apparent horizons and they could be naked. We also 
studied the future-directed outgoing null geodesies from 
singularities. Then we showed the existence of radial null 
geodesies from singularities in odd and even dimensions 
separately. This is a clear violation of cosmic censorship 
conjecture (at least the strong version). As we mentioned 
in the introduction, singularities cannot be naked in more 
than six-dimensional spacetime in the naive dimensional 
extension of the Einstein gravity. Therefore the Love- 
lock corrections worsen the situation in the aspect of the 
CCC. These results coincide with those of the spherical 
collapse of inhomogeneous dust in the Gauss-Bonnet the- 
ory, which is special case of the Lovelock theory, and the 
null dust collapse in the Lovelock theory. 

Compared with the Gauss-Bonnet gravity, we found 
that the Lovelock term changes the nature of singular- 
ities more drastically. In the Gauss-Bonnet case, only 
massless naked singularities could form in more than six 
dimensions. However, we found that massive naked sin- 
gularities can form in all odd dimensions in the Lovelock 
gravity. 

We also examined the strength of the singularities and 
found that i?pyfc^fc'^ , with the null tangent fc'^ , diverges 
at the singularities for all cases considered here. We also 
showed that the behavior of divergence depends on the 
spacetime dimensions in even dimension. However, it 
does not depend on the spacetime dimensions in the odd 
case. It is a characteristic feature of the Lovelock gravity. 

In this paper we had several assumptions which should 
be relaxed in future studies. First of all, we focused 
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on the marginally bound cases with W{r) = 1. But, 
our analysis can be easily extended to the nonmarginally 
bound cases (W(r) 7^ 1). We assumed that the mat- 
ter is a dust cloud. From the standpoint of the CCC, 
we should deal with more generic matter fields. We also 
assumed here spherical symmetry. Finally, most of our 
work is just the local analysis of singularity. The global 
visibility of singularity is still unclear. These are some of 
the remaining issues which will be hopefully addressed in 
the near future. 



On the other hand, using the metric of the outer region 
is also rewritten as 

d,' - - (^f{R^)tl dt' + Rldnl (A4) 

Of course, they should be identical and then we see that 

2 \ 



1 = f{R^)n 



2 Ry. 



f{R^) 



(A5) 
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Appendix A: Junction between inner and outer 
solutions 



In this appendix, we will address the junction between 
the inner and outer solutions. We suppose that the inner 
solution is the collapsing dust cloud presented in the main 
text and the outer solution is a static vacuum black hole. 
We follow the argument in Refs. [l^H^], where they dis- 
cussed the similar junctions in the Gauss-Bonnet theory 
and Einstein theory. 

Let S to be the boundary of the two regions. The inner 
solution to Eqs. (gOl) and is 



ds^ = -df + R'^dr^ + R^dnl, 
and the outside solution to Eq. (j4]) is 

1 



rfs^ ^ -f{R)dT^ 



fiR) 



dR^ + R^dni 



(Al) 



(A2) 



We suppose that S is described by parametric equations 
R = Rs(t) and T — T^{t). In addition, it is natural to 
think that the boundary is the comoving, that is, r = 
ro =constant because of dust. The induced metric on S 
from the metric of the inner region is written as 



holds. 

The extrinsic curvatures of S evaluated from inner 
metric are 



-Kl = 0, 



(A6) 
(A7) 



In terms of the metric of the outer region, the extrinsic 
curvatures also have another expression as 



• 2 \ -1 



(A8) 
(A9) 



where 



'f + R^ = fT. (AlO) 

The continuity of the extrinsic curvature implies us 

/3=1. (All) 

From the definition of /3, we have the following two equa- 
tions 

"^^'^ f{R^)-\ (A12) 



dt 



dR^ 
~dt 



where 

Then we see that 



(A13) 



(A14) 



(A15) 



is required. 

It is remembered that is the solution to Eq. ^ as 
in the form of 



ds' 



-df + R^dVll 



(A3) 



(A16) 
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and R is the solution to Eq. ((33 

k 



(\ 2m 



(A17) 



As seen easily, Eq. (IA15[) is satisfied only if 

/X = Firo). (A18) 

This shows that the inner region can be naturally joined 
to the outer region if = F{ro) is satisfied. 

Appendix B: Homogeneous dust collapse 

In this section, we briefiy discuss homogeneous case 
which is just a special case of inhomogeneous collapse. 
In this case, the spacetime metric becomes 



(Bl) 



Equation (|33)) reduces to 



k 

m— 1 



2m 



M 



(B2) 



where M is a constant. From this equation, we can see 
that the behavior of a near the singularity (a = 0) as 



a{t) ^{ts-t) 

ait) - t)"/"+i 



(in odd dimensions) (B3) 
(in even dimensions) (B4) 



where ts stands for the epoch of the singularity forma- 
tion. From these behaviors, it turns out that singularity 
is ingoing-null in odd dimensions and is spacelike in even 
dimensions according to Ref. 3l|. In the Gauss-Bonnet 
homogeneous case ^3|, singularity in five dimension is 
ingoing-null and is spacelike in dimensions higher than 
six. Therefore, as in the main text, we can confirm that 
the Lovelock terms changes the nature the singularity. 
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